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To perform a local orbital analysis of electronic and magnetic interactions, we construct the 
Wannier functions (WFs) of the Fe 3d orbitals in the parent compound of the recently discovered 
iron pnictide superconductors, LaFeAsO, and a comparison material LaFePO. Comparing the WFs 
for the stripe antiferromagnetic order with those for no magnetic order, the difference is a significant 
spreading ("rfelocalization") of specifically the d^y and dxz (but not dy^) WFs, where parallel Fe spins 
lie along the x direction. The WF basis gives a tight-binding representation of the first principles, 
density functional based Fe-derived bands. Comparing hopping parameters, it is found that changes 
due to stripe antiferromagnetism, even if it is weak, enables more isotropic hopping involving spin- 
majority electrons in the Fe 'idxz (but not the 3dyz) orbital. This change, counterintuitively, actually 
reinforces electronic anisotropy. Further insight is gained by comparing the WFs of LaFeAsO and 
LaFePO, identifying how the difference in WFs is related to the difference in hopping integrals 
and showing how the pnictide atom is influential in forming the stripe antiferromagnetism. Kinetic 
energy considerations suggest that orbital fluctuation, in addition to spin fluctuation, may contribute 
to the decrease in observed ordered moment compared to the calculated values. 



I. BACKGROUND AND MOTIVATION 

Since the first report from Hosono's group^ of su- 
perconductivity at Tc— 26 K in F-doped LaFeAsO, 
hundreds of experimental and theoretical papers 
on these iron-pnictide compounds have appeared, 
aimed at elucidating various properties, including 
synthesizing new compounds to achieve higher T^, 
measuring basic quantities (e.g. magnetic suscepti- 
bility, NMR, ARPES), and modeling and simulating 
to obtain explanations and predictions. Thanks to 
these efforts, there are now several families of these 
iron pnictide superconductors, including the 1111- 
family (e.g. LaFeAsO, CaFeAsF), 122 family (e.g. 
BaFe2As2), Ill-family (e.g. LiFeAs) and a more 
complicated 22426-family (e.g. Fe2As2Sr4Sc206), 
with Tc up to 56 K.[2| Several aspects have been 
clarified: the superconductivity lies in primarily 
iron 3d bands'^ and is not phonon-mediated;'* the 
ground state in most classes is a stripe antifer- 
romagnetic phase with a significantly reduced Fe 
magnetic moment compared to theoretically calcu- 
lated valuej^i^ it is a moderately correlated sys- 
tem where a Coulomb interaction U«3 eV might 
be appropriately There is discussion that the su- 
perconducting order parameter may have a new s± 
character. ^-^ 

Despite a great deal of progress in under- 
standing the electronic structurei2r— and magnetic 
interactions f^^— some basic questions remain unre- 
solved. One of them is: what is the underlying mech- 
anism of the structural transition from tetragonal 



to orthorhombic in the parent compounds of iron- 
based superconductors? This question is especially 
challenging in the 1111-compounds (e.g. LaFeAsO), 
where the structural transition is observed (as the 
temperature is lowered) to occur before^ the mag- 
netic transition (from nonmagnetic to stripe antifer- 
romagnetic order which we denote as Qm AFM). It 
would have been natural to think that the stripe an- 
tiferromagnetic ordering of Fe provides the driving 
force for the structural transition because it intro- 
duces electronic anisotropy. (Table III in reference 
[17[ provides a summary of the structural transition 
temperature Tg and stripe antiferromagnetic tran- 
sition temperature T^r of several iron pnictide com- 
pounds.) 

Noting that the structural transition and mag- 
netic transition occurs simultaneously in the 122 
compounds (e.g. BaFe2As2), a possible argument 
is that the magnetism is in fact present, in the 
form of medium-range order, antiphase boundaries, 
etc., near the structural transition but its detection 
is greatly suppressed by strong spatial or temporal 
fluctuation. The suggestion by Mazin and Johannes 
that magnetic antiphase boundaries may be the 
dominant cxcitatiorii^ has already stimulated nu- 
merical estimations by the present authors.^- With 
a time resolution of 10""'^^ s, photoemission experi- 
ments by Bondino et al.^^ implied a dynamic mag- 
netic moment of Fe with magnitude of 1 fis in the 
nonmagnetic phase of CeFeAsOo.soFo.n, which is 
comparable to the ordered magnetic moment of Fe 
in the undoped antiferromagnetic CeFeAsO com- 



pound. The fluctuation strength should be much 
stronger in 1 Ill-compounds than 122 compounds 
based on the fact that the measured Fe ordered 
magnetic moment in 1111-compounds (^ 0.4 (Ib) 
is much less than in 122 compounds (^ 0.9 fis) and 
they are much smaller than DFT predicted value 
(~ 2 fig)^^ One factor is that interlayer coupling 
of FeAs layers is stronger in 122 compounds than 
1111-compounds because the interlayer distance in 
122 compounds (~6 A) is significantly smaller than 
1111-compounds (~8-9 A),H Interlayer interaction 
should help to stabilize the ordered Fe magnetic mo- 
ment by reducing fluctuations (reducing two dimen- 
sionality). 

In this paper we address the effect of magnetic 
order, and of the pnictide atom, on the strength, 
character, and spin-dependence of Fe-Fe hopping 
processes by using a Wannier function represen- 
tation based on all five Fe 3d orbitals, and only 
these orbitals. Several previous studies of the elec- 
tronic structure have pointed out some aspects of 
the influence of the pnictide, or chalcogenide, atoms 
(due to size or chemical identities) and also of their 
positionsj^iSi"— We provide one example of the effect 
of the pnictogen atom (comparing LaFaAsO with 
LaFePO) in this paper, where the effect of the pnic- 
togen is included precisely but indirectly through the 
Wannierization process. This allows us to present 
results in an Fe-centric picture. This local orbital 
representation provides insight into both electronic 
and magnetic behavior even when the fundamental 
behavior is primarily itinerant. 
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II. CALCULATIONAL METHODS 

We begin with first principle calculations us- 
ing the full-potential local-orbital code^ (FPL08) 
with local density approximation (LDA) exchange- 
correlation (XC) functional^ (PW92) and the ex- 
perimental lattice constants and internal atomic co- 
ordinates for the compounds LaFeAsO and LaFePO, 
as given in our previous work^ i^^i^^ . To obtain a con- 
sistent local orbital representation and the resulting 
hopping amplitudes, we then construct real-space 
Wannier functions (WFs) derived from Fe 3d orbitals 
in both NM and Qm AFM phases. The WFs used in 
this paper, as implemented in the FPL08 code, are 
constructed by projecting the Bloch functions from 
a specified energy range onto chosen atomic orbitals, 
roughly following the method of Ku et ali^^^ The 
resulting Wannier orbitals retain a symmetry that 
is common to both the atomic orbital and the point 



FIG. 1: LaFeAsO band structure with highlighted Fe 
Sdyz and 'idxz fatband characters in the NM (top panel) 
and Qj\/ AFM (bottom panel) phases. Compared to the 
NM phase, the Fe id^z bands near Fermi level in the Qm 
AFM phase, especially along T — X and V — Y directions, 
change dramatically due to the formation of the stripe 
antiferromagnetism with large ordered Fe magnetic mo- 
ment of 1.9 IJbB- 



group symmetry of the site. These WFs provide an 
explicit basis set of local orbitals that give a tight 
binding representation, complete with on-site ener- 
gies and hopping amplitudes to neighbors as distant 
as necessary to represent the chosen bands. In this 
paper we project onto the conventional real Fe 3d 
orbitals, with the energy range corresponding to the 
region with strong Fe 3d character in the bands. 
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FIG. 2: LaFePO band structure with highlighted Fe 
3dyz and Sd^z fatband characters in the NM (top panel) 
and Qm AFM (bottom panel) phase. Compared to 
LaFeAsO, the Fe Sd^z bands near Fermi level in the Qm 
AFM phase change less significantly from the NM phase, 
due to the relatively small ordered Fe magnetic moment 
of 0.5 hb- 



III. DIFFERENCES IN BAND 
STRUCTURES 



The differences in electronic structure that we will 
emphasize result from the changes due to stripe mag- 
netic order, and the differences between LaFeAsO 
with larger ordered moment, and LaFePO, with 
smaller calculated moment (experimentally nonmag- 
netic). The necessary band structures are shown in 
Fig. [1] for LaFeAsO and H for LaFePO, where in 
each case the Fe 3dyz and Sd^z characters are high- 
lighted. The total energy of LaFePO, which is exper- 
imentally found to be nonmagnetic, is only slightly 



lower (2 meV/Fe) in the Qi\i AFM phase than the 
nonmagnetic phasej^ii so the incorrect prediction for 
LaFePO is actually a fine detail, and suggests it is 
nearly antiferromagnetic. The calculated Fe mag- 
netic moment of LaFePO in the Qm AFM phase 
is 0.52 hb- In LaFeAsO, the calculated moment is 
near 1.9 fj,B, substantially larger than the measured 
value of 0.36 /is as has been widely discussed (see, 
for example, Refs. [slillj). 

For our calculations and discussion we have chosen 
the a;-axis along the direction of aligned Fe spins, as 
shown in Fig. [Sj the corresponding zone boundaries 
are denoted X and Y in the band plots. The nonmag- 
netic band structures of the two compounds are very 
similar, differing only in some fine details that do not 
arise in our analysis. However, the band structures 
in the Qm AFM phase of the two compounds differ 
substantially, which is due to the difference in the 
Fe magnetic moment (1.9 vs. 0.5 IjLb)^'^- The sim- 
ilarities and differences provide a way to study the 
effect of magnetic order in these compounds, and 
specifically to show that even small magnetic order 
has substantial consequences. Since the AFM and 
NM phases in LaFePO are nearly degenerate, our 
results have relevance to the effect of (longitudinal) 
magnetic fluctuations of the Fe atom. 

The panels in Fig. [T] illustrate that the magnetic 
order substantially simplifies the band structure very 
near the Fermi level, which is all near F in this dou- 
bled (magnetic) cell. The other difference to notice 
is the great difference in band structure along F-X 
and F-Y directions. Figure [5] shows the influence of a 
weak stripe antiferromagnetism (0.5 (Ib) on the non- 
magnetic band structure. The overall band structure 
remains the same except for some bands near the 
Fermi energy, where the main change is the separat- 
ing of the Fe Sdxz bands away from the Fermi level, 
which causes disappearance and change of topology 
of certain pieces of the Fermi surface of the Fe 3dxz 
bands. Note that the Fe 3dyz bands change insignif- 
icantly, leaving the bands near the Fermi level dom- 
inated by Fe Sdyz character. This difference indi- 
cates that even a weak stripe antiferromagnetism 
has a very strong symmetry breaking effect on the 
3dxz and 3dyz bands, which are equivalent in the 
nonmagnetic state. As a result, even a weak stripe 
antiferromagnetism induces a large anisotropy, let 
alone the much stronger (calculated) antiferromag- 
netism in FeAs-based compounds. (The much bigger 
anisotropy in the stripe AFM phase in LaFeAsO is 
evident by comparing Fig. [T]and[2J) 
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is often called "orbital ordering," but based on the 
calculated populations, discussed below, this is more 
properly thought of as an itinerant cousin of orbital 
ordering.) txy denotes the hopping parameter of the 
dxz—dxz hopping in the y direction, and tyx the dyz — 
dyz hopping in the x direction. In the nonmagnetic 
case, by symmetry 



'^xy — '^yx — ^5 
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whereas they differ in the Qa/ AFM state, t^x will 
denote dxz — dxz hopping in the x direction, and sim- 
ilarly tyy denotes dyz— dyz hopping in the y direction 
(see Fig. [I. 

Let U and U' denote the intra-orbital and inter- 
orbital Coulomb repulsion, and Jh the inter-orbital 
Hund's exchange constants. Our purpose is to es- 
timate the difference in kinetic energy gain of the 
two configurations shown in Fig. [31 At the level of 
second-order perturbation theory, the kinetic energy 
dyz hopping in the x direction 
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T gain from the dyz 

spin down (Fig- 13^) is 



FIG. 3: (color online) Possible orbital orderings of iron in 
iron-pnictides. Left panel: Both (a) and (b) form the Qm 
AFM ordering. However, (a) is favored because it gains 
more kinetic energy from nearest-neighbor hoppings ac- 
cording to second-order perturbation theory (see text). 
Right panel (from top to bottom) shows the simplified 
symbols for Fe 3dyz and Sd^z orbitals, the chosen x and 
y directions, up arrows for spin up electrons and down 
arrows for spin down electrons, where black arrows for 
'idyz orbital and red arrows for Sd^z orbital. 



^Eyx = -tlx/iU' - Jh). 



A similar kinetic gain of 



^Exy — 



<y/{U' - Jh) 



(2) 



(3) 



dxz hopping in the y direction 
are much smaller and can be 



comes from the dxz 

(Fig. HJi). txx and tyy 

neglected (see Table |l|. Therefore, the total energy 

gain from NN hopping of Fig. [3^ is 



AE{a) = AExy + AEyx - -2tV(C/' - Jh), (4) 
while it is 



AE{b) ^ -2t^/U 



(5) 



IV. 3d,z 

REPOPULATION 



Due to the strong influence of stripe antiferro- 
magnetism on the band structure (even when weak 
as in LaFePO), the orbital distinction and repop- 
ulation of the Fe idxz and idyz electrons suggests 
various means of analysis. The strong intra-atomic 
anisotropy discussed above is sometimes referred to 
as orbital ordering, but with the orbital occupations 
far from integers, the anisotropy also has a substan- 
tial itinerant (collective) component. Here we con- 
sider briefly the alternative, local viewpoint. 

Figure |3] shows two underlying (idealized) orbital 
populations, both of which are consistent with the 
Qa/ AFM symmetry. (This orbital differentiation 



for Fig. [3b. Because U is larger than U' — Jh, 
the orbital ordering in Fig. ^i&) is favored over Fig. 
[3ljb) by kinetic fluctuations. This result is a more 
transparent form of an analysis presented by Lee et 



V. TIGHT BINDING HOPPING 

PARAMETERS AND WANNIER 

FUNCTIONS 

Figure [4| shows the Wannier functions (WFs) of 
all five Fe 3d orbitals in both NM and Qa/ AFM 
(majority spin) phases of LaFeAsO, using the same 
value of isosurface in all cases. In the NM (spin- 
degenerate, tetragonal) phase, all five WFs for Fe 



TABLE I: The hopping parameters (in eV) of the Pel 
nearest neighbor Fe2 and Fe4 atoms and next-nearest 
of LaFeAsO. The highlighted (itahcized and boldface) 



neighbor Fe3 atom in the nonmagnetic and Qm AFM phases 
entries are discussed in the text. 
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3d orbitals have their density strongly concentrated 
on the Fe site. All Fe minority spin 3d WFs in the 
Qa/ AFM phase remain almost the same as in the 
NM phase, so they are not shown. The majority 
spin WFs for 3dyz, 3dx^^y2 and 3(i^2 orbitals remain 
very similar to the corresponding Wannier functions 
in the NM phase, as can be seen in Fig. |4l The 
significant change is that the majority spin WFs for 
3dxz and 3dxy become more delocalized in the Qm 
AFM phase, with significant density at the nearest 
neighbor As sites, the effect being especially large for 
the dxz orbital. This difference reveals that the ma- 
jority spin Fe 3dxz and 3dxy orbitals mix much more 
strongly with nearest-neighbor As Ap orbitals in the 
Qa/ phase than in the NM phase. The AFM or- 
der involves a highly anisotropic magnetization, and 
resulting difference in majority and minority poten- 
tials, that produces this strongly orbital-dependent 
effect. 

Using these WFs as the basis gives a tight bind- 
ing representation, for which the hopping parame- 
ters are obtained from matrix elements of the Wan- 
nier Hamiltonian using the FPL08 code. The corre- 
sponding band structures of LaFeAsO and LaFePO 
are already shown in Fig. [T] and Fig. [5] and the re- 
sulting tight binding bands (not shown) fit very well 
the corresponding DFT-LSDA Fe-derived bands in 
both NM and stripe AFM phases. 

TableUpresents the hopping parameters of the Fel 
3dyz, 3dxz, and 3dxy orbitals to all the 3d orbitals 
of its nearest neighbor Fe2 and Fe4 atoms and next 



nearest neighbor Fe3 atom in LaFeAsO. (See Fig. [3] 
for the definition of each Fe atom.) 

The on-site energies (in eV) of all the five 3d or- 
bitals in the NM phase and Qa/ AFM phase in both 
LaFeAsO and LaFePO are shown in Table |lll In 
the Qa^ AFM phase, the on-site energies are shown 
separately for both spin up (majority spin) and 
spin down (minority spin) orbitals. For LaFeAsO, 
the dxz , dyz energies lie at the Fermi level, ^^2 and 
dx2-y2 lie 0.1-0.3 eV below while dxy lies about 0.2 
eV above, without magnetic order. With AFM or- 
der, the minority on-site energies do not change 
greatly, whereas the majority levels fall by 0.7-0.8 
eV, thereby affecting the hybridization with the As 
Ap orbitals. The changes in LaFePO are smaller, 
corresponding to the smaller (factor of ^3) magnetic 
moments. 

The hopping parameters reported here are similar 
to the corresponding hopping parameters reported 
by Lee et al^ and Haule et al^. (The differences 
between our results and those of Lee et al. reflect 
the fact that, although the original bands are the 
same and the Wannier transformation is formally 
the same, the Wannier transformation is not unique 
and depends somewhat on some details in the im- 
plementation.) Our hopping amplitudes are not di- 
rectly comparable to those reported by Cao et al^ 
who focused on the hopping parameters from As Ap 
orbitals to Fe 3d orbitals and to its nearest neigh- 
bor As Ap orbitals. As shown in Table HI in the 
NM phase, txy — tyx » txx — tyy, which indicates 
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TABLE II: The on-site energies (in eV) of the 4^2, 
d^2_y2, dyz, dxz, and d^y Fe orbitals in the NM and 
Qa/ AFM phases in LaFeAsO and LaFePO. In the Qa/ 
AFM phase, the on-site energies are shown separately 
for the spin up (majority spin) and spin down (minority 
spin) orbitals. 
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FIG. 4: LaFeAsO Wannier functions of Fe 3d orbitals in 
the NM phase (top panel), and those for the majority 
spin in Qm AFM phase (bottom panel). The Wannier 
functions of Fe 3d orbitals for the minority spin in the 
Qa/ AFM phase remain almost the same as in the NM 
phase. The important difference to be observed is that in 
the Qm AFM phase, the (majority spin) Wannier func- 
tions of 3da;z and Sd^y orbitals (and only these) are more 
extended, with much increased density at neighboring As 
sites. The isosurface has the same value (density) in each 
panel. 



that the hopping (through As atoms) of dxz — dxz 
{dyz— dyz) in the y (x) direction of the electrons in Fe 
idxz {idyz) orbital is favored over the x (y) direction. 
The hopping process for Fe Zdxz (Sdyz) electrons is 
anisotropic. Global tetragonal symmetry is retained 
because the Fe 3dxz and Sdyz electrons hop in dif- 
ferent directions, which enforces the equivalence of 
the X and y directions. 

In the Qm AFM phase, the corresponding hop- 



ping parameters (both spin up and spin down) are 
either the same or very close to the NM value, ex- 
cept for two cases. These differences are intimately 
related to the changes in the corresponding WFs, as 
we now explain. The first one is the dxz — dxz hop- 
ping between parallel spin Fe neighbors {x direction) 
of a majority spin electron, whose absolute value in- 
creases significantly from the NM case (from -0.06 
to -0.29 eV, see the highlighted numbers in Table IJ). 
This opens an extra hopping channel in addition to 
the original dxz ~ dxz hopping in the y direction. In 
the NM state, the electrons in the dxz or dy^ orbitals 
separately only hop in one direction (in the sense 
that the hopping parameters in other directions are 
relatively small) . The dramatic change of the 3dxz 
bands near Fermi level from NM to Qm AFM, noted 
in several previous studies, can be traced to this dif- 
ference. 

The other case is for dxy — dxy hopping, again 
between parallel spin atoms {x direction). In the NM 



phase, the dxy 



hoppings in x and y directions 



are the same by symmetry, with an amplitude of 
0.18 eV. In the Qa/ AFM phase, this hopping in the 
y direction for both spins is slightly enhanced to 0.23 
eV. However, the dxy— dxy hopping in the x direction 
is significantly enhanced to 0.33 for the majority spin 
and suppressed to 0.07 for the minority spin. These 
differences shows that the broken symmetry has a 
strong effect on the dxy orbital's environment. 

The magnitude of the changes of the hopping pa- 
rameters in the two special cases mentioned above, 
and thus the magnetic order induced changes in 
WFs, is directly related to the magnitude of the or- 
dered Fe magnetic moment in the Qm AFM state, 
which is evident by comparing the case of LaFeAsO 
and LaFePO (see Table HI and HIT)) . The iron atom 
in the Qm AFM state in the former compound has 
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TABLE III: The hopping parameters (in eV) of the Fel 3dyz, Sd^z, and Sd^z orbitals to all the five 3d orbitals of its 
nearest neighbor Fe2 and Fe4 atoms and next-nearest neighbor Fe3 atom in the nonmagnetic and Qm AFM phases 
of LaFePO. 
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a large ordered magnetic moment of 1.9 /is while 
in the latter compound it is very weak, only 0.5 
fjLB^ in DFT-LSDA calculations. The difference in 
the ordered Fe magnetic moment is consistent with 
the change of hopping parameters of d^z — dxz and 
dxy — dxy in the x direction of the spin majority elec- 
tron from the NM to the Qj\/ AFM state, as shown 
in Table Hand Table |llll 

The difference in the changes of the hopping pa- 
rameters of each Fe Sd orbital from NM phase to 
Qm AFM phase is related to the spin polarization 
of each orbital in the Qm AFM phase, as shown in 
Table IIVI The id^z orbital has the largest moment 
(0.51 hb in LaFeAsO), followed by the 2>dxy orbital 
(0.48 fiB in LaFeAsO). The other three orbitals have 
significantly smaller moments (less than 0.41 /is in 
LaFeAsO). It is clear that the orbital with larger or- 
bital spin magnetic moment has bigger changes in 
the relevant hopping parameters. The difference of 
the relevant hopping parameters between LaFeAsO 
and LaFePO can also be traced to the difference in 
the orbital spin magnetic moment. 

The transition to the Qa/ AFM state is accompa- 
nied, in a local picture and to second order, by an 
extra kinetic energy gain of 



l^Ex 



<x/iU' - Jh) 



(6) 



from the hopping process of dxz — dxz hopping in the 
X direction, which is comparable with AExy. (Note 
that AExx is negligible in the NM state.) A substan- 
tial extra kinetic energy gain can also be obtained 



from the dxy — dxy hopping in the x direction. The 
anisotropy arises because the majority-spin electron 
in the 3dxz orbital can hop in both directions {i.e. to 
both parallel and antiparallel spin neighbors), while 
others in the Sdxz and 3dyz orbitals can basically 
only hop in one direction. This anisotropy is re- 
flected in a large symmetry lowering of the Sdxy or- 
bital in the AFM phase. The anisotropy leads to 
a large spin polarization (orbital spin magnetic mo- 
ment) in the Sdxz and Sdxy orbital, which may also 
be related to the tetragonal to orthorhombic struc- 
tural transition such that the lattice constant along 
the aligned-spin direction (x direction in this paper) 
becomes shorter than the other direction {y direction 
in this paper, thus a < b). 



The additional Sdxz — Sdxz hopping and the en- 
hancement of the 3dxy — 'idxy hoppings, both in the 
X direction of the spin majority electron, promote 
kinetic energy gain. However, as pictured in Fig. 
[3^, the 3dxz spin up electron of Fel atom cannot 
hop in the x direction due to the Pauli principle. 
In order to take advantage of this extra kinetic en- 
ergy gain of AExx, the spin up occupation number 
of the 3dxz orbital should not be unity but instead 
must fluctuate. The same situation happens to the 
Sdxy orbital. The competition between the kinetic 
energy gain and Pauli principle results in a reduced 
magnetic moment and is possibly one mechanism of 
orbital fluctuation. 



TABLE IV: Occupation numbers and spin polarizations in 3d orbitals in the NM and Qm AFM phases of LaFeAsO 
and LaFePO compounds. Sn is the difference of the total occupation number in each orbital between the Qm AFM 
phase and the NM phase, m is the spin magnetic moment in each orbital in the Qm AFM phase. 





LaFeAsO 


LaFePO 




NM Qm 

up dn 5n m 


NM Qm 

up dn 5n m 


z' 

2 2 

X -y 
yz 

xz 
xy 


0.71 0.89 0.48 -0.05 0.41 
0.57 0.80 0.45 0.10 0.34 
0.65 0.85 0.57 0.11 0.28 
0.65 0.86 0.35 -0.10 0.51 
0.68 0.88 0.39 -0.11 0.48 


0.69 0.75 0.64 0.01 0.11 
0.54 0.59 0.50 0.01 0.09 
0.67 0.72 0.64 0.01 0.08 
0.67 0.75 0.56 -0.03 0.19 
0.67 0.71 0.62 -0.00 0.09 



VI. SUMMARY 

In this paper, we have compared the electronic 
structures of LaFeAsO and LaFePO, in both NM 
and Qj\/ AFM phases, and find that the stripe an- 
tiferromagnetism affects very differently the various 
Fe 3d orbital characters, even when the stripe an- 
tiferromagnetism is weak. By comparing LaFeAsO 
to LaFePO (and looking at similar results for other 
1111 and 112 compounds^), we find that the pnic- 
tide atom and the structure are influential in the 
formation of Qm AFM phase, consistent with sev- 
eral earlier reports that did not provide any de- 
tailed analysis. This information was obtained from 
a tight-binding representation for Fe Sd electrons 
based on first principles Wannier functions. 

In the nonmagnetic phase the electrons in Fe Sd^z 
and idyz orbitals have very different amplitudes to 
hop in the x and y directions, resulting from the po- 
sitions and chemical character of the pnictide atoms. 
Anti-intuitively, this "anisotropy" is almost gone for 
majority spin electrons in the AFM phase, when the 
3dxz (or 3dyz) electron can hop equally to parallel 
and antiparallel neighbors (both x and y directions) . 
This change is accompanied by a lowering of sym- 
metry, and extension in space, in the Sd^y Wannier 
function. The (large) changes in the near neighbor 



hopping parameters of the Sd^z and Sd^t/ orbitals 
in the x direction is directly connected to the much 
larger orbital spin magnetic moments of these two 
orbitals than the other three orbitals. 

The anisotropy in hopping in the Fe 3dyz, Sd^^j, 
and 3dxy orbitals also favors orbital fluctuation by 
providing extra kinetic processes, which are partly 
compensated by the Pauli principle which inhibits 
the hopping processes, and which we expect to en- 
hance fiuctuations in the corresponding orbital occu- 
pation numbers (orbital fluctuation). Such fluctua- 
tions would reduce the ordered Fe magnetic moment 
in the Qm phase, bringing them closer to the ob- 
served ordered moments. On the other hand, inter- 
layer hoppings of the Fe 3d electrons in the z direc- 
tion may help to stabilize the Fe magnetic moment 
in the Qm AFM phased 
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VIII. CALCULATIONAL METHODS 

We perform first principle calculations using 
the full-potential local-orbital code^^ (FPL08) 
with local density approximation (LDA) exchange- 
correlation (XC) functionaP'5 (PW92). In the cal- 
culations, we use experimental lattice constants 
and internal atomic coordinates for the compounds 
LaFeAsO, LaFePO, CaFeAsF, SrFeAsF, BaFe2As2, 
SrFe2As2, and CaFe2As2, as used in our previ- 
ous work^iiiii^. For the other two hypothetical 
compounds LaFeNO and LaFeSbO, the lattice con- 
stants and internal atomic coordinates are taken 
from the optimized equilibrium values of first prin- 
ciple calculations'^^ done in the Qm AFM phase us- 
ing GGA (PBE) XC functional^, since such calcu- 
lations were proven to predict the correct equilib- 
rium lattice constants and internal atomic coordi- 
nates compared to the experimental values in all the 
known iron pnictide compounds] ^^'^^ 



IX. FURTHER OBSERVATIONS 

Similar hopping parameters compared to 
LaFeAsO have been obtained for the CaFeAsF, Sr- 
FeAsF, and MFe2As2 (M=Ba, Sr, Ca) compounds, 
(which have similar FeAs layers), as shown in Table 
IVl However, replacing As in LaFeAsO with other 
pnictides (N, P and Sb) results in similar t^y, tyx 
and tyy but different t^x- Compared to LaFeAsO, 
the txx for the majority spin electron in the Q7\/ 
AFM phase is reduced for LaFeNO and LaFePO, 
but enhanced in LaFeSbO. The importance of the 
pnictide for the formation of the Qm AFM phase is 
evident. 

Another important factor is the interlayer hop- 
pings. The interlayer distance of FeAs layers in 1111- 
compounds is in the range of 8.2 -9.0 A and it is 
much smaller in 122-compounds, ranging from 5.9 
A to 6.5 A. The interlayer hopping parameters of 
Fe 3d electrons in the z direction are negligible in 
1111-compounds but become substantial for certain 
hoppings in 122-compounds, especially in CaFe2As2, 
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TABLE V: The ratios of hopping parameters ta;y, tyx, 
txx and tyy in the NM and Qm AFM phases of a few 
iron-pnictides. 







yz 


xz 


compound 
(mag. mom.) 




NM Qa/ 

up dn 


NM Qm 

up dn 


LaFeNO 
(1.86 iib) 


^yx/^xy 
^yy I'^xx 


-0.30 -0.33 -0.27 
-0.03 -0.05 -0.05 


-0.30 -0.31 -0.31 
-0.03 -0.14 0.06 


LaFePO 

(0.52 hb) 


^yx/y^xy 
^yy I'^xx 


-0.37 -0.37 -0.34 
-0.09 -0.09 -0.09 


-0.36 -0.36 -0.36 
-0.09 -0.15 -0.03 


LaFeAsO 
(1.90 iib) 


^yx/^xy 

tyy ftxx 


-0.33 -0.42 -0.29 
-0.06 -0.09 -0.09 


-0.33 -0.35 -0.35 
-0.06 -0.29 0.09 


LaFeSbO 

(2.45 iib) 


tyx/txy 
tyy /hxx 


-0.26 -0.39 -0.21 
-0.07 -0.11 -0.11 


-0.26 -0.28 -0.27 
-0.07 -0.38 0.16 


CaFeAsF 
(1.75 iib) 


tyx/txy 
tyy/txx 


-0.36 -0.43 -0.34 
-0.06 -0.08 -0.08 


-0.36 -0.37 -0.37 
-0.06 -0.27 0.08 


SrFeAsF 
(1.96 iib) 


tyx/txy 
tyy /I'xx 


-0.35 -0.43 -0.31 
-0.08 -0.10 -0.10 


-0.35 -0.37 -0.37 
-0.08 -0.31 0.08 


BaFe2As2 
(1.88 iib) 


tyx/txy 

tyy l%xx 


-0.32 -0.40 -0.29 
-0.08 -0.10 -0.10 


-0.32 -0.34 -0.34 
-0.08 -0.28 0.07 


SrFe2As2 
(1.78 hb) 


tyx/txy 
tyy ftxx 


-0.33 -0.40 -0.31 
-0.08 -0.10 -0.10 


-0.33 -0.34 -0.35 
-0.08 -0.28 0.06 


CaFe2As2 
(1.67 ilb) 


tyx/txy 
tyy jtxx 


-0.33 -0.38 -0.32 
-0.08 -0.10 -0.10 


-0.33 -0.35 -0.35 
-0.08 -0.28 0.06 



whose interlayer distance of FeAs layers is only 5.9 A. 
Certain interlayer hopping parameters are as large 
as 0.15 eV for idxy and 'id^^ orbitals, and 0.07 eV 
for idyz, 3dxz and Sdx^-y^ orbitals, calculated in the 
Qm AFM phase for CaFe2 As2 , which has the small- 
est interlayer distance. 

The large interlayer hopping parameters for the 
Fe 3dxy orbital, which at first sight seems very sur- 
prising, becomes clear by noting that the 3dxy Wan- 
nier orbital is strongly distorted from its symmetric 
atomic shape to its nearest neighbor As atoms above 
and below the Fe plane, as shown in Fig. 4 in the 
original paper. This extension in the z direction will 
favor interlayer hoppings, especially when the inter- 
layer distance is small, as in the case of CaFe2As2. 
For comparison, the interlayer hopping parameters 
(if not zero) are less than 0.01 eV in LaFeAsO. The 
increasing hopping of Fe 3d electrons in the z di- 
rection increases the interlayer coupling, and may 
inhibit fluctuations and thereby help to stabilize the 
ordered Fe magnetic moment in the Qm AFM phase. 
The kz dispersion correlates with the experimen- 
tal observations that the measured Fe magnetic mo- 
ments in the Qm AFM phase are significantly larger 
in 122-compounds {^ 0.9 hb) than 1111-compounds 
(-0.4/xb). 



